Journal of Biomolecular Structure and Dynamics

ObjectiveThe complement regulatory protein CD59 has been shown to promote SNARE complex assembly, yet its interaction with vesicle-associated membrane protein 2 (VAMP2) remains poorly characterized. This study aims to identify the key domain in VAMP2 that mediates the CD59 interaction and to evaluate whether CD59 point mutations affect SNARE complex assembly. MethodsThe interaction between CD59 and VAMP2 was examined by immunofluorescence confocal microscopy and co-immunoprecipitation (co-IP). The effect of CD59 on SNARE complex assembly was assessed by co-expressing CD59 with the three core SNARE proteins (syntaxin-1, SNAP-25, and VAMP2) and detecting complex formation by western blotting. Four CD59 single-point mutants were generated and evaluated in SNARE assembly assays. AlphaFold3 was employed to predict the interaction between CD59 and individual VAMP2 domains (confidence threshold: ipTM + pTM [≥] 0.75). Truncated VAMP2 constructs were further characterized by molecular dynamics simulations and co-IP. Results(1) CD59 directly bound VAMP2 and promoted SNARE complex assembly without altering individual SNARE protein levels. (2) All four CD59 single-point mutants retained the ability to promote SNARE assembly at a level comparable to wild-type CD59, despite showing differential effects on binding stability in molecular dynamics simulations. (3) The proline-rich (P-rich) N-terminal domain of VAMP2 was identified as the key binding interface; its deletion abolished the CD59 interaction, whereas deletion of the SNARE motif did not. ConclusionCD59 promotes SNARE complex assembly through interaction with the P-rich N-terminal domain of VAMP2. The examined point mutations do not impair this function, suggesting that these sites may tolerate substitutions or that redundant contact residues maintain the interaction.

Cancer is a serious public health problem, and is becoming more common, with a projected increase in deaths and more than 25 million new cases in coming decades. A number of molecular mechanisms are involved in the tumoral process, with one of them, never in mitosis A-related kinase 2 (NEK2), a serine/threonine protein kinase, being a frequent target of amplification in various malignancies that is responsible for chromosomal instability, aneuploidy and activation of several oncogenic pathways. Available kinase inhibitors are not yet optimized in terms of their pharmacokinetic properties for clinical use, and current therapies, such as chemotherapeutic agents or immunotherapies are often limited by their resistance. In silico methods represent an effective tool to search for novel potent inhibitors, before testing in animals, with time constraints and limited resources. To find new inhibitors of NEK2, we used E pharmacophore-based modeling and structure based virtual screening in this study. NEK2 was chosen as the target for therapeutic intervention and an energy optimized pharmacophore model was employed to screen the Enamine REAL library of millions of compounds. Pharmacodynamic and Pharmacokinetic properties of the Top hits were tested using ADMET profiling. These were further screened using molecular docking (standard precision and extra precision) and virtual screening to obtain three lead compounds 1, 2, and 3 which have docking score of -7.414, -8.037 and -7.562 respectively. MM-GBSA calculations were used to estimate the binding free energies for these complexes, which were determined to be -54.92, -54.18 and -49.23 kcal/mol. Lastly, 100 ns molecular dynamics simulations have been run to evaluate complex stability in dynamic situations. The overall results of the MD showed the overall stability of the NEK2-ligand complexes, and thus these three compounds are promising NEK2 inhibitor candidates and could be further validated in vitro and in vivo for clinical application. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=133 SRC="FIGDIR/small/742111v1_ufig1.gif" ALT="Figure 1"> View larger version (70K): org.highwire.dtl.DTLVardef@14e9cfeorg.highwire.dtl.DTLVardef@24ec78org.highwire.dtl.DTLVardef@20be75org.highwire.dtl.DTLVardef@1b80c02_HPS_FORMAT_FIGEXP M_FIG C_FIG

Engulfment and Cell Motility protein 1 (ELMO1) regulates cell migration, phagocytosis, and cytoskeletal remodeling, positioning it as a compelling therapeutic target across kidney diseases, oncology, enteric infections and inflammation. Despite this potential, no approved therapeutics or clinically validated small-molecule modulators of ELMO1 currently exist. ELMO1 functions by forming a complex with DOCK180 (or DOCK2) to activate the small GTPase Rac1, and the recent structural resolution of the ELMO1/DOCK2 complex now provides an opportunity to target this protein-protein interface directly. Here, we present the first investigation into the druggability of the ELMO1/DOCK2 complex and report the initial virtual screening to identify small-molecule inhibitors of this interaction. Molecular dynamics (MD) and free energy level (FEL) studies were carried out to validate the potential of the predicted hits. This work establishes a computational framework for the development of the first generation of ELMO1-targeted therapeutics. In addition to demonstrating the drugability of ELMO1, this work introduces two open-source Python tools for the rapid analysis and visualization of protein-protein interaction and ligand-protein MD trajectories from DESMOND output files. These tools are designed to be broadly accessible, offering practical utility to the wider DESMOND user community.

Nipah is a virus borne zoonotic disease that mostly occurs in Southeast Asian countries. Nipah virus (NiV) infection is potentially lethal as it can cause fatal encephalitis, and as of now, there are no vaccines or therapeutic treatments for this disease. The highly pathogenic NiV belongs to the Paramyxovirus family and the natural host of NiV is fruit bats of the Pteropodidae family. Though NiV has not been detected in the US, rather recently the Camp Hill virus (CHV) in shrew species was identified in Camp Hill, Alabama. This virus is characteristically similar to Nipah, and belongs to the henipavirus genus that include the Nipah and Hendra viruses. NiV has two surface glycoproteins namely the attachment glycoprotein (NiVG) and the fusion glycoprotein (NiVF). Being associated with host cell attachment and fusion of viral membrane, these glycoproteins tend to facilitate the viral entry to the host cell. The present paper describes the structures of the aforesaid viral glycoproteins based on the available database of protein structures and examines if they have an evolutionary link to other viral glycoproteins. Certain similar glycoproteins found in different countries are also explored. The structure-function relationship for these two NiV glycoproteins are determined using time-based molecular dynamics (MD) simulation. The interfacial ligand-receptor interactions between NiV proteins with human receptor and monoclonal antibody, Fab 66, are examined. A primary motivation for these structural studies is that the results might potentially contribute to future therapeutic developments and drug designs to effectively combat NiV.

Blue light using Flavin (BLUF) proteins are microbial photoreceptors that are involved in various physiological responses. Their occurrence and biochemical properties in fungi remain poorly understood. Here, we investigated a putative BLUF photoreceptor from the corn-smut fungus Mycosarcoma maydis (MmBLUF). Domain analysis, multiple sequence alignment of BLUF core regions, and structural modelling indicated conserved canonical BLUF fold and flavin-pocket residues. However, when heterologously expressed, UV-visible and fluorescence spectroscopy revealed different spectral behaviour than canonical BLUF protein. Further, we tested the role of extended N-terminus in modulation of chromophore binding by expressing N-terminus truncated protein variants. Our results suggest that the unusual spectral behaviour is not linked to the truncation construct (extended N-terminus), which also showed similar spectral features, indicating that the extended N-terminus is unlikely to account for an unusual photodynamics characteristics. Our findings support MmBLUF as a structurally conserved putative fungal BLUF-like photoreceptor with different photochemical properties. Further studies are required to establish its chromophore identity, photocycle and function of this unusual BLUF-like domain from fungal system.

Group A Streptococcus (GAS), or Streptococcus pyogenes is almost exclusive and greatly adapted human pathogen. It causes a wide array of clinical symptoms, ranging from minor infections of the skin and soft tissues to pharyngitis, meningitis, pneumonia, bacteraemia, cellulitis, puerperal sepsis, and necrotising fasciitis. The risk of S. pyogenes infection is known to be influenced by several host characteristics, including age, underlying diseases like diabetes, varicella, or skin lesions, both chronic and acute, and certain risk behaviours such as use of drugs. Household size and overcrowding are two environmental factors that significantly affect the transmission of S. pyogenes. The majority of cases occur spontaneously in the community, and preventative opportunities are still limited. A large portion of GAS-related mortality is found in low-income areas and communities. Based on aforementioned public health risk, the creation of effective therapeutic vaccines would be an excellent addition to current control measures. The purpose of this work is to address the need for new instruments to aid in the elimination of S. pyogenes infections. The discovery of high antigenic regions in several highly conserved proteins brings us one step closer to developing peptide vaccines capable of influencing the different phases of S. pyogenes infection, providing more effective defence and greater serotype coverage. This study used various techniques of immunoinformatics to design an effective multi-epitope vaccine that produced neutralising antibodies against multiple strains of S. pyogenes.

The binding affinities of a chimeric analog of a riboflavin derivative linked to biotin, (6- (7,8-dimethyl-2,4-dioxo-3,4-dihydrobenzo[g]pteridin-10(2H)-yl)hexyl 5-((3aS,4S,6aR)-2- oxohexahydro-1H-thieno[3,4-d]imidazol-4-yl)pentanoate), referred to as C6-Rf-biotin-tag, to the riboflavin binding retain or streptavidin are in the M range, 1.29 {+/-} 0.277 and 3.00 {+/-} 0.459, respectively. These values suggest that C6-Rf-biotin-tag has potential applications in diagnostic assay and labelling target flavin binding proteins. The C6-Rf-biotin-tag which was characterized with respect to physical and biochemical properties retains UV/Vis spectroscopic and fluorescence behavior similar to riboflavin.

Autophagic clearance of hyperphosphorylated tau is impaired in tauopathies, leading to the progressive accumulation of toxic tau species. Fungal metabolites provide a rich yet largely untapped source of neuroactive molecules with therapeutic potential. Here, we investigated metabolites from Lions Mane (Hericium erinaceus), Magic Mushrooms (Psilocybe spp.), and Ergot fungi (Claviceps spp.) using a computational drug-discovery workflow. We characterised their structural diversity, predicted blood-brain barrier permeability and toxicological properties, and integrated network pharmacology with protein-protein interaction and functional enrichment analyses to identify autophagy-related targets. Peroxisome proliferator-activated receptor gamma (PPARG), glycogen synthase kinase-3 beta (GSK3B), and casein kinase 2 alpha 1 (CSNK2A1) emerged as the three most promising candidates, given their complementary roles linking autophagy and tau pathology and their attractiveness as targets for multi-target drug discovery. Their interactions with fungal metabolites were evaluated by molecular docking, Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) rescoring, molecular dynamics simulations, complemented by machine-learning quantitative structure-activity relationship (QSAR) modelling as an additional, ligand-based line of evidence. Molecular dynamics and MM/GBSA analyses confirmed stable, target-specific binding for Corallocin A and Erinacerin M (PPARG), Chaetopyranin and Ergocryptine (GSK3B), and Hericioic Acid D and Isohericerin (CSNK2A1), alongside Emodin, a reference compound with previously reported activity against all three targets. QSAR predictions were informative primarily for the PPARG candidates, which fell within the models applicability domain; predictions for the remaining candidates fell outside their respective models applicability domains and were therefore not interpretable as evidence for or against their prioritisation. Reanalysis of an independent hippocampal proteomic dataset from Alzheimers disease patients showed CSNK2A1 protein levels to be significantly altered in the CA3 subfield, providing an additional, correlative line of support for this target; PPARG and GSK3B showed no significant changes in protein abundance, which does not preclude their functional involvement given their extensive post-translational regulation. Overall, these findings identify fungal metabolites as promising multi-target autophagy modulator candidates and provide a systematic computational strategy for prioritising them for experimental validation in tauopathies.

Valproic acid/sodium valproate (VPA) is a widely prescribed anticonvulsant and has also been used against certain tumor cells. It is a potent modulator of gene expression. Its ability to induce apoptosis has been well documented in HeLa cells. However, another form of cell death - mitotic catastrophe - has not yet been explored in VPA-treated HeLa cells. Here, we investigated the effects of VPA treatment on mitotic catastrophe characteristics, including morphological features and their frequencies, fluorescence intensity signals of caspase-2 and p53, and the expression and abundance of DNMT1 and DNMT3B. An increased frequency of mitotic catastrophe was observed not only morphologically, but also through enhanced induction of caspase-2, involvement of p53, at least under more drastic VPA treatment, but without a decrease in DNMT1 or DNMT3B levels. Additionally, enhancement of mitotic catastrophe coincided with a reduction in mitotic chromosome abnormalities. Increased DNMT3B expression following VPA action, may be favored by previously reported chromatin decondensation induced by this drug. Enhanced CpG methylation of specific DNA sites could thus be promoted. In conclusion, VPA was shown to trigger metabolic pathways linked to different forms of cell death in HeLa cells, supporting its oncosuppressive potential.

The Dengue virus (DENV), a major global health concern, causes dengue fever, predominantly affecting tropical and subtropical regions. The NS2B-NS3 protease complex of DENV is critical for viral replication, making it a promising target for antiviral drug development. This study investigates the structural dynamics of the NS2B-NS3 protease under varying pH conditions, integrating experimental and computational approaches. The recombinant NS2B-NS3 protease was expressed in E. coli, purified using affinity chromatography, and analysed for purity via SDS-PAGE. Circular dichroism (CD) spectroscopy revealed pH-dependent secondary structural changes, indicating stability at neutral to slightly basic pH and destabilization under acidic and highly basic conditions. Dynamic light scattering (DLS) analysis demonstrated protein aggregation and structural heterogeneity under extreme pH levels. Complementary in silico techniques, including homology modelling and molecular dynamics (MD) simulations, provided detailed insights into the conformational changes of the protease. The modelled structure, validated and refined through computational tools, revealed structural stability at physiological pH, with notable disruptions at pH extremes. DSSP (Dictionary of Secondary Structure of Proteins) and principal component analyses highlighted significant secondary structural transitions, especially at acidic pH where -helices and {beta}-sheets transformed into random coils Docking and MD simulation was carried out with the small molecule for therapeutic analysis between Dengue NS2bNS3 and small molecule. This study emphasizes the pH-dependent conformational plasticity of the NS2B-NS3 protease, contributing to the understanding of its functional mechanisms and providing a foundation for the rational design of pH-specific inhibitors. These findings underscore the importance of structural biology in advancing therapeutic strategies against dengue fever.

Carbohydrates are essential molecules involved in a wide range of biological processes. While their structural representation in the Protein Data Bank (PDB) has been hindered due to their highly diverse and flexible nature, advancements on X-ray crystallography and CryoEM are accelerating the obtaining of their experimental structures, which in turn supports improved analyses of their structural properties, such as the role of their non-chair conformations on biological processes. In this context, the current work provides a comprehensive and up to date overview of carbohydrate structures within the PDB through classification, quantification, and characterization of monosaccharides and glycan chains, incorporating Cremer-Pople puckering coordinates to evaluate pyranose conformations. Our findings show that, even with a steep increase in CryoEM structures including more complex oligosaccharides, carbohydrate data remains concentrate dinsmallglycanchains, with a about 4% of non-chair conformation sintopresolution structures. Free energy landscapes for the pyranose ring in the presence and absence of nine different binding proteins, including both enzymes and non-enzymes, reinforce the preferred stability of chairs and the rarity of conformational induction upon binding.

The spread and evolution of antimicrobial resistance (AMR) pose significant threats to public health, food safety, and sustainable development. Metallo bea lactamases (MBLs), a class of carbapenemases that hydrolyze carbapenems, often considered last resort antibiotics, are of particular concern and underscore the urgent need to develop effective inhibitors. Recent advancements in bicyclic boronate inhibitors represent significant progress toward this goal. Although the mechanistic basis of NDM1 inhibition by bicyclic boronate inhibitors remains incompletely elucidated, this study systematically characterizes the binding profiles and intermolecular interaction landscapes between NDM1 and four cyclic boronic acid based inhibitors, vaborbactam (RPX7009), taniborbactam (VNRX 5133), xeruborbactam (QPX7728), and ledaborbactam (VNRX-5236). We established the NDM-1-mediated hydrolysis of two representative carbapenem antibiotics, meropenem and imipenem, as the reference system for comparative analysis. While conventional docking workflows often struggle to model metalloenzyme active sites accurately, the AMDock (Assisted Molecular Docking) platform, integrated with AutoDock optimized zinc force field, enables robust and reliable docking simulations for this enzyme family. Our AMDock based molecular docking not only identified the critical NDM1 binding site residues governing hydrogen-bond networks with ligands, but also revealed that bicyclic boronate inhibitors show substantially greater binding stability than both carbapenems and their monocyclic boronate counterparts. A follow-up 100 ns all atom molecular dynamics simulation in GROMACS provided atomistic level evidence for this enhanced performance: RMSD, RMSF, radius of gyration, continuous hydrogen bonding occupancy, and 2D free energy landscape analyses all consistently supported the higher binding affinity of bicyclic boronates. Ultimately, this study characterizes two promising drug candidates, QPX7728 and VNRX5133, as competitive inhibitors that can effectively rescue carbapenem activity against NDM1 MBL mediated antimicrobial resistance. This visualization elucidates interactions between NDM1 and its two ligand classes, carbapenem antibiotics and cyclic boronate inhibitors. It reveals the stability and interaction strength underlying the distinct inhibition pattern of bicyclic boronate inhibitors in the enzyme active site, offering a promising strategy to combat antibiotic resistance.

Interleukin-13 (IL-13) is an immunomodulatory cell signaling cytokine that has been implicated in neurodegenerative disease and chronic inflammation. IL-13 binds to its low and high affinity receptors, IL-13 receptor 1 (IL-13R1) and IL-13 receptor 2 (IL-13R2), respectively, with residence times that vary accordingly. As the binding kinetics of the cytokine-receptor complexes influence cellular responses, we employed the molecular dynamics (MD) simulation-based{tau} -random acceleration molecular dynamics method ({tau}RAMD) to compute relative residence times for wild-type (WT) IL-13 and 19 IL-13 mutants to the two receptors. Comparison with experimental kinetic data shows that the{tau} RAMD computations capture the trends in residence times. Analysis of simulated dissociation trajectories of the cytokine-receptor complexes reveals two distinct dissociation pathways of IL-13 from each of the receptors. This study thus pinpoints key determinants of the interaction of IL-13 with its receptors which could be targeted for therapeutic design. Statement of SignificanceCytokines are regulatory proteins that bind to cell surface receptors and thereby send signals to the cellular interior. Interleukin-13 (IL-13) is a cytokine that has a low and a high affinity receptor. It has important physiological roles, and its deregulation is involved in diseases such as atopic dermatitis and asthma. We employed a molecular dynamics simulation-based method to compute the effects of changes in the sequence of IL-13 on the lifetimes of complexes of IL-13 and its receptors. Comparison with experiments supports the validity of the computational approach and analysis of the simulations reveals two distinct ways in which IL-13 dissociates from each receptor. These results thus provide a map for targeting IL-13 - receptor interactions for the design of therapeutics.

Drug targets and vaccines for bacterial diseases are challenging to develop. Computational subtractive genomics and reverse vaccinology offer a viable substitute for conventional drug and vaccine development approaches. In recent years, Pasteurella multocida, a Gram-negative pathogen responsible for severe infections in humans along with major diseases in livestock and poultry, has developed resistance to multiple drugs. As a result, there is an urgent demand for research into reliable, efficient, highly potent, and cost-effective strategies to combat the rise of Pasteurella multocida infections. To address this, the present study integrated an in silico subtractive genomics and reverse vaccinology workflow to pinpoint proteins with potential value as both drug and vaccine targets. From the complete Pasteurella multocida Pm70 proteome, we identified 446 proteins that had no similarities with human proteins, out of which 129 were found to be essential for the pathogens survival. These 129 proteins are therefore potential candidates for drug development. Further metabolic pathway analysis, virulence potential, subcellular localization, and structural analyses revealed two outer membrane proteins, namely, 47 kDa outer membrane protein and lipopolysaccharide-assembly protein LptD, as promising candidates for therapeutic development. Subsequent antigenicity analysis revealed both outer membrane proteins as potential vaccine targets, exhibiting highly favorable antigenicity scores. Candidate B- and T-cell epitopes were systematically mapped from these two targets based on antigenicity, immunogenicity, toxicity, and allergenicity potential and only high-confidence epitopes were assembled using adjuvants and class-specific linkers to construct the final multi-epitope vaccine (MEV). The designed MEV construct was demonstrated to be highly antigenic, non-toxic and non-allergenic to humans, while possessing favorable structural and physicochemical properties, and exhibited strong binding interactions with human toll-like receptors with minimal molecular deformability. Although further experimental validation is recommended to fully assess their efficacy, the present study establishes a solid foundation for the development of novel therapeutics and vaccine candidates against P. multocida, ultimately contributing to the reduction of its rising infections.

In TMPRSS2 fusion-positive prostate cancer, ERR is involved in regulation of ERG and promotes the androgen receptor independent signaling in the cancer progression. The ERR binds to the ERREs (estrogen-related receptor response elements) present at -5042 bp of the TMPRSS2- promoter and enhances the ERG overexpression that causes prostate cancer progression. To dissect the structural basis of the ERR recognition to the TMPRSS2 promoter DNA, we have purified the full-length ERR (ERRFL), NTD deleted construct (ERR{Delta}NTD), and the DNA-binding domain (ERRDBD) proteins and performed the binding analysis with 30 bp TMPRSS2-promoter DNA (5' -AGTCCAAGGTCGGTGGATC ACAAGGTCAGG-3'). Circular dichroism analysis showed that all three ERR proteins adopt native secondary structures. DNA binding induced subtle changes in the secondary structures, while enhancing the thermal stability (Tm) of all ERRa proteins. Binding analysis showed that ERRDBD bound weakly to the DNA, whereas ERRFL and ERR{Delta}NTD exhibited substantially higher affinities ~120-fold and ~131-fold than ERRaDBD, respectively. Small-angle X-ray scattering (SAXS) analyses revealed a dimeric ERRFL structure and an ERRFL-DNA complex (2:1) structure in solution and fitted well with Alpha Fold model of apo and DNA bound complex of ERRFL. Furthermore, 100 ns dynamics simulations on apo and DNA-bound ERRa proteins showed that all proteins remained structurally stable, with flexibility largely confined to loop regions of ERRa proteins. Our biophysical, DNA binding and structural analyses have revealed the mechanism involved in ERR recognition of the TMPRSS2- promoter DNA, which provides insight into ERR-mediated transcriptional regulation and development of anticancer drugs against ERR-driven prostate cancer.

Cytosolic phospholipase A2 (cPLA2), encoded by PLA2G4A plays a crucial role in platelet activation through thromboxane A2 biosynthesis. Defects in this pathway can lead to bleeding disorders; however, the precise role of cPLA2 in Glanzmann thrombasthenia (GT) or GT-like syndromes remains poorly understood. The current study emerged from a fortuitous discovery of novel PLA2G4A variants, p.Asp455Val (D455V) and p.Glu119Lys (E119K), among 3 GT patients from two unrelated families in south India. We envisaged a hypothesis to identify a hitherto unknown role of cPLA2 that could be attributed to GT-like bleeding symptoms in these patients. An extensive computational workflow, including variant classification, conservation analysis, stability prediction, molecular docking, molecular dynamics simulations, principal component analysis, and binding free energy (MMPBSA) calculations, were employed to assess mutation-induced alterations in the cPLA2 protein. The current findings depicted that both variants affected conformational dynamics, ligand recognition patterns, and local interaction networks, without causing any drastic global structural destabilization. In particular, the D455V partially disrupted the catalytic domain, potentially reducing substrate (PAPC) stabilization and hampering interactions with catalytic products. Likewise, the E119K variant partially affected the C2 domain, potentially altered electrostatic interactions and enhanced regional flexibility, which may compromise product (AA) stabilization. Further analyses revealed mutation-induced changes in collective motions, protein-ligand interactions, solvent accessibility, and secondary structural elements. Overall, our findings promote the hypothesis that both E119K and D455V variants may partially alter phospholipid-associated interactions, and cPLA2-mediated downstream signaling. Moreover, this study underscores cPLA2 as a potential modifier of platelet function and bleeding severity.

Oncogenic mutations are known to occur at non-random rates and specific hotspots, yet the structural and thermodynamic factors underlying these hotspots remain poorly understood. This study presents molecular dynamics simulations of 22 cancer hotspots across 10 oncogene families, simulated as histone-DNA complexes. Across 18 of 22 mutations, thermodynamic effects were predominantly localized to within 6 [A] of the DNA-histone interface (mean capture 104%), with van der Waals interactions driving the effect at the thermodynamic extremes, indicating that oncogenic mutations alter nucleosome binding through precise local contact changes rather than global structural rearrangements. Further, analysis revealed a gene-family correlated pattern in thermodynamic stability. RAS family mutations showed a consistent trend toward nucleosome stabilization (mean {Delta}{Delta}G = -4.82 kcal mol-1), while kinase domain mutations trended toward destabilization (mean {Delta}{Delta}G = +40.12 kcal mol-1), a difference reaching statistical significance in this exploratory analysis (Mann-Whitney U, p = 0.008). This interface-specific mechanism, combined with the gene-family-correlated thermodynamic pattern, provides a biophysical framework for understanding nucleosome-level contributions to cancer hotspot mutation biology.

Per-and polyfluoroalkyl substances (PFAS) constitute a diverse class of persistent synthetic chemicals utilized across industrial, medical, and consumer sectors that are pervasive global pollutants. Exposure to PFAS is linked to adverse impacts on both innate and adaptive immune systems. Human lactoferrin (hLF) is a key antimicrobial component of the developing innate immune system present in colostrum and breast milk. We hypothesized that hLF is a potential PFAS binding protein related to PFAS immunotoxicity. The results of thermal stability experiments indicated that all 11 tested PFAS bind and destabilize the structure of hLF. Notably PFBA, PFOS, HFPO-DA, and 6:2 FTSA decreased apo-hLF melting temperatures from 64oC to [≤] 37oC, suggesting that PFAS exposures destabilize the native hLF protein under physiological conditions. Relative binding affinities (Kd) ranged from 0.2-11 mM across tested PFAS. Molecular docking was used to confirm experimental binding affinities and identify molecular interactions involved with PFAS binding. Calculated Gibbs Free Energies of binding ranged from -4.4 to -8.8 kcal/mol. Together, these results demonstrate that PFAS bind hLF at affinities comparable to human serum albumin and other PFAS binding proteins, and that some PFAS can destabilize hLF protein structure at physiologically relevant temperatures and conditions.

It is essential to experimentally evaluate antiviral efficacy predicted by molecular docking and related in silico approaches, which this study achieves for three bisphosphonates (alendronate, minodronate, zoledronate) and treprostinil, identified as potential antiviral therapies in our previous work. These investigational compounds had no detectable cytotoxicity at concentrations up to 25M in Vero E6 cells but also failed to show any antiviral effect against the PQ.8.1 isolate of SARS-CoV-2, when compared to control drugs (ensitrelvir, nirmatrelvir, and remdesivir). It appears that for bisphosphonates, any protective effect is likely due to other virus- or host-mediated mechanisms of action. This study demonstrates the importance of iterating between theoretical hypotheses and experiments to elucidate underlying mechanisms at play.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease characterised by progressive motor neuron degeneration. Mutations in the SOD1 gene represent the second most common genetic cause of ALS (ALS), and distinct SOD1 missense variants present with markedly different clinical profiles. A4V leads to an aggressive form of the disease (median survival [~]1y), H46R confers a mild, slowly progressive course and I113T exhibits an intermediate phenotype. The molecular basis by which these mutations produce divergent clinical outcomes remains poorly understood. We performed extensive classical molecular dynamics simulations of wild-type SOD1 and the three ALS-associated variants in the apo monomeric state to attempt to investigate the mechanisms behind such phenotypic differences. Structural stability, global compactness, and conformational flexibility, as well as analysis of collective motions between residues and estimation of free energy, were assessed. The H46R, A4V, and I113T variants exhibited distinct dynamic behaviours, highlighting differences in structural stability, local flexibility, and intramolecular interactions. These findings suggest that specific structural regions may contribute differently to protein dysfunction and could represent key elements for understanding the relationship between molecular dynamic properties and the differing clinical severity associated with these variants. Most strikingly, H46R exhibited exceptional structural stability across every analytical level, the lowest global deviation, most attenuated local flexibility, strongest internal dynamic coordination, and the deepest, most confined free energy basins of any system examined. This convergent multi-layered evidence of structural restraint provides a compelling mechanistic basis for the mild and slowly progressive clinical course of H46R ALS, suggesting that enhanced conformational rigidity, rather than bulk destabilisation, is the defining biophysical feature of this variant, and that its pathogenic mechanism operates through a route fundamentally decoupled from the aggregation-driven toxicity that characterises the more aggressive SOD1-ALS mutations.